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O Abstract 

■ 

^3 ■ We compute the photon spectrum and the rate for the decays B{D) — > lv\^ 

These photonic modes constitute a potentially large background for the purely 
Qh. leptonic decays which are used to extract the heavy meson decay constants. 

<-j \ While the rate for D — ► luj is small, the radiative decay in the B meson 

case could be of comparable magnitude or even larger than B — > [iv. This 
J^j ' would affect the determination of fs if the r channel cannot be identified. We 

j_j ■ obtain theoretical estimates for the photonic rates and disscuss their possible 

experimental implications. 
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I. INTRODUCTION 



The leptonic decays of heavy mesons are of great interest both theoretically and ex- 
perimentally. The purely leptonic decay B~ — > V~V\ can be used to extract \V U \\ by using 
predicted values for the B meson decay constant fs from lattice calculations. Or conversely, 
if \V u b\ is measured in semileptonic decays, the leptonic decay is in principle the only mode 
to access On the other hand, the situation in the charmed mesons is clearer given our 
knowledge of the CKM mixing angles involved (|V^| and |V^|). Recently the observation 
of the decay D s — > /iz/ M has been reported [1J. Although the results are still preliminary 
it is a first step towards understanding the behavior of heavy meson decay constants. The 
experimental difficulty in the measurement of purely leptonic decays of heavy pseudoscalars 
is due mostly to the well known effect of helicity suppression: back-to-back leptons must 
make a spin final state, but the anti-neutrino is right-handed and forces the charged lepton 
to this helicity which introduces a factor of the lepton mass in the amplitude. In the end 
the decay rate is suppressed by the factor (m;/m#) 2 , where ran is the heavy meson mass. 
For instance for the B meson 

nB - m = ^Wfs (—) 2 ml (l - 4) (1) 
87T \m B J \ m B J 

where the pseudoscalar meson decay constant is defined by 

(Olgj^l -7 5 )b\B{P)) = if B P li (2) 

and analogous expressions can be written for D mesons. Thus, when the charged lepton is 
an electron the purely leptonic decay is practically inaccessible. At the other extreme, when 
the charged lepton is a tau there is no suppression. However the observation of this decay 
is experimentally difficult (only one visible particle in the final state). Muons seem more 
promising to allow the observation of the decay constants of D and B mesons. However in 
B — > \w the suppression factor is still 4 x 10~ 4 giving a branching fraction of 1CT 7 — 1CT 6 . 

There are other decays that indirectly involve the heavy meson decay constants. For 
instance, the decay B — > -nlvi is expected to be largely dominated by a B* pole diagram at 
very low recoiling pion energies |2|-44|]. This implies the presence of the vector- meson decay 
constant /g» which can be related to Jb by Heavy Quark Spin Symmetry |§. However this 
region of phase space is difficult to access due to kinematic suppression. 

In this letter we investigate the decay modes B~(D~) — ► l~ui"f. There are two types of 
contributions: Internal Bremsstrahlung (IB) and Structure Dependent (SD) photon emission 
||. As is known the IB contributions are still helicity suppressed. On the other hand, the 
SD contributions are reduced by the electromagnetic coupling constant a but they are not 
suppressed by the charged lepton mass. Therefore, what in principle could be regarded as a 
mere radiative correction to the purely leptonic decays has the potential to be of comparable 
magnitude and in some cases even much larger. In what follows we analyze the situation in 
charmed and beauty mesons. In Sec. 2 we establish the phenomenological relevance of these 
decays and in Sec. 3 we discuss theoretical estimates of the unknown constants involved. 
Conclusions and comments are presented in Sec. 4. 
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II. PHENOMENOLOGY OF H — > /ivy 

As mentioned in the previous section, the [i modes are the most interesting from the 
point of view of the extraction of heavy meson decay constants. We will concentrate on the 
case I = but the treatment for / = e is analogous and the numerical differences between 
these two cases will be stressed when relevant. 

The emission of a real photon in leptonic decays of heavy mesons can proceed via the 
two mechanisms mentioned in Sec. |I[ The possible IB diagrams are shown in Fig. 1. The 
corresponding amplitude is given by 

M PB = mJ B x { (i^L _ Ml _ J5)u 

where P, pi and k are the four-momenta of the decaying meson, the charged lepton and the 
emitted photon, respectively, e is the polarization of the photon and F is the electromagnetic 
form-factor of the B meson. (This is not simply related to the Isgur-Wise function J7| given 
that the matrix element of the electromagnetic current between two B mesons receives 
contributions from both the heavy and light quarks). The important feature of (f|) is that it 
is suppressed by a factor of the lepton mass, as one would expect of bremsstrahlung photons. 

On the other hand, the SD diagrams in Fig. 2 involve the contributions from heavy 
intermediate states coupling to the initial heavy pseudoscalar and the photon. In Fig. 2 
we show the contributions from vector and axial-vector mesons. The helicity suppression is 
avoided because the meson directly coupling to the lepton pair has spin one. These types of 
diagrams were previously considered in the context of light pseudoscalar decays, in particular 
in 7r — > Ivy and K — > Ivy M. It was also shown there that this phenomeno logical picture in 
fact correctly accounts for the SD contributions. 

The general form of the amplitude corresponding to Fig. 2(a) requires the knowledge of 
the B*B , y coupling. This is defined in the process B* — > Bj with the amplitude given by 

A(B* -> Bi) = enve^e^vpkjXg (4) 

where v is the velocity four- vector of the decaying particle, k is the photon four-momentum, 
A is the polarization four-vector of the B* and \xy is a constant characterizing the strength 
of the Ml transition. Eqn. (|]) can be used to write the amplitude of Fig. 2(a) by having 
first the B meson decaying into a real photon and an off-shell vector meson, which after 
propagating decays weakly through 

(0\qj li (l-7s)b\B*(v 1 X))=if B .€ fl (5) 

Eqn. (|5|) defines the vector-meson decay constant. Inserting the heavy vector-meson propa- 
gator i(—g^ v + v fl v u )/2(v.k + A), where A = m B * — m B , we obtain its contribution to the 
SD diagram 

A V (B - iivi) = efMv —^— fB^elv^Cs (6) 
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(3) 



with Cs = Uijs(l — 75)1*1/ the lepton current. 

There will also be contributions from heavy axial-vector meson states (J p = 1+). The 
Heavy Quark Symmetries M of the strong interactions and their consequences in the hadron 
spectrum will help us to identify the relevant contributions. In the Heavy Quark Limit the 
spin of the heavy quark decouples from the light degrees of freedom. Thus their total angular 
momentum j\ can be used as a good label. The meson angular momentum is J = 1/2 ±j;. 
For the ground state ji = si, the spin carried by the light degrees of freedom and therefore 
a (0~, l - ) spin doublet is predicted. This can be identified with the (B, B*). If we allow for 
orbital angular momentum L — 1, then there will be two additional spin doublets as excited 
states. They correspond to ji = 1/2 and ji = 3/2 and their parity is even: (0 + , 1 + ) and 
(l + ,2 + ). Thus in principle there could be two different axial-vector mesons contributing to 
B — > fiuj. We write down the amplitude for the process in Fig. 2(b), which is the generic 
contribution of an axial-vector meson intermediate state. The coupling is of the form 

Af{B l -> Bj) = en\{v ■ kX ■ e — v • eX ■ k) (7) 

where the superscript i = 1/2, 3/2 identifies the spin-parity doublet to which the axial- vector 
meson belongs. The axial-vector decay constants f 1 / 2 and f^( 2 are defined analogously to 
Jb in (|5|). Putting these together the contributions of Fig. 2(b) take the form 

A<£ (B - yvj) = 2{v 6 ^ A Ai) [vek,-v ke,] C» (8) 

where Ai/ 2 = msj — tjib and A 3 / 2 = mB 1 — are the excitation energies of the ji = 1/2 
and ji = 3/2 even parity doublets and we again neglected pieces proportional to the charged 
lepton mass. It is worth noticing that the axial-vector pseudoscalar mass differences A1/2 
and A3/2 are quantities that remain constant in the limit tjiq — > 00 whereas A goes to zero 
in the same limit. In the Heavy Quark Effective Theory (HQET) flO| a term that goes as 
A 2 can be neglected in the B* propagator on the basis of the heavy mass suppression. In 
the case of the B 1 or B[ propagators in Fig. 2(b) the quadratic terms in Ay 2 and A 3 / 2 can 
still be considered numerically small even when they are quantities of order 1 in the heavy 
masses, when compared with 2m b- 

Most importantly, Heavy Quark Spin Symmetry can be used to obtain relations between 
the decay constants of particles in the same parity multiplet. For instance 

f B , = m B f B (9) 

and analogous relations can be written between the decay constants of the members of the 
jf = l/2 + and the jf = 3/2 + doublets. However HQS does not predict relations among 
decay constants of members of different spin-parity doublets. Then the contributions from 
excited states will imply the presence of new unknowns other than the pseudoscalar decay 
constant. Defining x = 2E 1 /rriB and y = 2Ei/rriB as the rescaled photon and charged lepton 
energies in the B rest frame the double differential decay rate is given by 

dxdy- 32^ Qm ^ X + + VTx + ^ 



x 



mg 1 \ l rag 

y\l -x)+ y(3x - x 2 - 2) + i(3s 2 - 4x - x 3 + 2) 



(10) 
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where we defined 



,.(*) fi 

7* = ^ (11) 
^VjB* 

as the relative axial-vector to vector meson coulping strength. The fact that there are no 
relations among decay constants of states of different spin-parity doublets will be reflected 
in the persistence of the unknown 7j's when we normalize to T(B — > [lis). We will analyze 
the 7j's and fiy from the theoretical point of view in the next section. For now let us assume 
that \iy and /J^ do not depend on the photon energy so we can integrate and compare with 
the purely leptonic decay. The result is 



DM 




6r a »v<[^f) x / dxx*{l-x){ - : 2A - + 1£ 1 }> (12) 



To have an idea of the potential importance of the photonic mode we take a definite value 
for the mass differences Aj = 600 MeV as suggested by the charmed meson system || . This 
gives approximately 

R B « 2/4 (l + (7l/2 + 73/2)2 j GeV 2 (13) 

which shows that unless there are unnaturally small B photon couplings (see next section) 
the photonic decay, which in principle could have been considered a small radiative correc- 
tion, will dominate the leptonic decay or will be at least of comparable magnitude. Therefore 
it is important to have a good theoretical understanding of the couplings of all the relevant 
intermediate states in order to substract these events as a background for B — > \w. The 
fact that, to this order in (m M /m^) 2 , the result of Eqn. (|T0|) is independent of the lepton 
mass implies that 

T(B -> evi) = T(B -> (jtisy) (14) 

and that 



T(B -> evi) > T(B -> ei/) (15) 

allows for the separation of both effects by using the electronic modes as well. Thus inte- 
grating (|H]) over the photon energy and integrating the resulting \x spectrum around the 
end point over a region the size of the experimental fi energy resolution in B — > \w will 
eliminate the background. Perhaps even more interesting, given the richness of the physics 
that enters in them, is the possibility of observing the photonic modes at branching ratios 
that will soon be accesible (~ 10~ 5 — 10 -6 ). 

The treatement of charmed mesons is entirely analogous to B mesons. With the obvious 
replacements in (|T2"D we obtain 
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ii^4x KT 2 /4 (l + (7l/2 + 2 73/2)2 ) GeV 2 (16) 

which is suppressed by the factor {m D /m B ) A in rescaling flT2|) . Therefore the effect is, as 
expected, less spectacular in the D mesons although of a branching fraction comparable to 
that of the effect in the B system. 

In the following section we discuss various theoretical estimates of \iy and the 7i's. 



III. THEORETICAL ESTIMATES 



The radiative leptonic decays depend crucially on the vector and axial-vector couplings 
to the heavy pseudoscalar, /iy and /i^ as well as the ratio of the axial-vector to vector 
meson decay constants. The photon couples to both the heavy and light quark pieces of 
the electromagnetic current. Heavy Quark Symmetry (HQS) fixes the heavy quark contri- 
bution. In the Heavy Quark Effective Theory |Hj] the coupling of the heavy quark Q to the 
electromagnetic field responsible for the B* to B transition is given by the operator 

QQ h v a^h v F^ (17) 



where h v is the heavy quark field charaterized by its four-velocity v and Qq is its electro- 
magnetic charge. This leads to 

1*™ = ^ (18) 

m Q 

where the superscript h indicates a contribution from the heavy quark. 

The piece of the B*B^ coming from the coupling to the light degrees of freedom is not 
given by HQS. In the SU(3) limit can be written as 1TB 



^ = Q q (3 (19) 

where Q q is the charge of the light degrees of freedom and (3 is an unknown quantity. In the 
non-relativistic quark model this is predicted to be 

Pnrqm = l/m q (20) 

with m q the constituent light quark mass (m u w uid ~ 330 MeV and m s w 450 MeV). On 
more general grounds (3 is expected to be of the order 1/Aqcd- If we take the HQS and 
the NRQM predictions for //W and fjfi respectively in order to have an estimate of the 
magnitude of the coupling we have the simple expression 

Qq + Q 1 (21) 

m Q m q 

which is the naive quark model result but now supplemented by HQS. In the charmed 
mesons, it is sufficient to explain the ratio T(D*° — > D 0r y) /T(D* + — > D +r y) [ ll|. Therefore 
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we will rely on fl2T|) for estimates in the B meson decays. For m& = 5 GeV and taking 
(3 = 3GeV _1 , we see that the value for the coupling 



-2GeV 



-i 



(22) 



is such that the ratio ( O ) could be sizeable. 

The axial- vector meson couplings are expected to be of the same order as \xy. In fact the 
NRQM gives = \iy j\P& and i^pj^ — ^/2/3/iy, where the superscripts indicate an axial- 
vector meson belonging to the (0 + , 1 + ) and (1 + , 2 + ) doublets repectively. The contributions 
of Fig. 2(b) are also proportional to f\, the axial-vector meson decay constants. But in the 
NRQM these are zero given that the wave function at the origin for an orbitally excited 
state vanishes due to the presence of a centrifugal barrier. Although this seems to imply 
that 7j <C 1, relativistic effects could be important and modify this prediction drastically. 
Unfortunately there is no conclusive experimental indication from other hadronic systems 
about the size of the axial-vector meson contributions. In experiments involving 7r — > fiw-f 
and K — > fiw-f the ratio of the axial- vector to vector contributions is consistent with zero 
and also with being large. 

In order to address corrections to this simple picture it is useful to consider an effective 
theory coupling heavy hadrons (in this case mesons) to goldstone bosons and also low energy 
photons. This theory incorporates both Heavy Quark Symmetry and Chiral Symmetry by 
introducing the spin-parity doublet in the form of 4 x 4 matrices as follows 
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^75} 



(23a) 



S n 



to - Bl) 



(23b) 



a 



(1 + 



1 



(23c) 



where a = 1,2,3 is the SU(3) index of the light degrees of freedom in the hadron. H a 
corresponds to the (0~,1~) ground state doublet; S a is the excited state doublet (0 + , 1 + ) 
corresponding to the spin of the light degrees of freedom being ji = 1/2 and T£ is the excited 
state doublet (l + ,2 + ) corresponding to ji = 3/2. 

An important set of corrections are those arising from the loop diagrams of Fig. 3, where 
the photon couples to a goldstone boson in the loop. This type of correction therefore is not 
suppressed by the heavy mass. They were first considered in Ref. Jl3| for the D*D^ coupling 
in the context of Heavy Hadron Chiral Perturbation Theory (HHChPTh) and introduce a 
non-analytic dependence on the quark masses, m l J 2 as well as on the mass differences A 
and Aj that break SU(3) in (|IDD. To calculate these diagrams we need the coupling of two 
heavy hadrons to one goldstone boson. These are given by |L2[ 



C ln = gTr H a H b fi ba 'j 5 + g'Tr SaS^ajs + g"Tr 
+f'Tr ISaTfA^] + f'Tr \H a S^ bal , 



(24) 
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with the traces taken over Dirac indices and where the first three terms correspond to the 
axial couplings between members of the same spin-parity doublet and the last two terms give 



the transitions between doublets. As it was pointed out in [121 the axial couplings between 



the (1 + , 2 + ) doublet and the ground state vanish to this order in the chiral expansion, that is 
Tr H a Tfr Afj.ba'-fn vanishes. This means that the diagrams in Fig. 3 corresponding to having 
the (Bi, B£) inside the loop can be neglected. On the other hand we can also see from (pl[) 
that to this order there will be no B[B7r coupling which prevents a contribution from the 
axial-vector meson to the loop. There are four diagrams remaining which correspond to the 
heavy mesons in Fig. 3(a) being (from left to right) (B,B*,B*); (B,B*,B[)\ (B,B%,B[) 
and (B, Bq, Bi). In terms of the couplings, the first one corresponds to a correction to 
the two following to a correction to y$ coming from B[ — ► B^y and the last one to a 
correction to /i^ 2 ' ) from B\ — > B*y. Finally the diagrams of Fig. 3(b)-(c) vanish in the limit 
mi — > and will be strongly suppressed. Therefore the only diagrams contributing are the 
ones represented in Fig. 3(a), with pion and kaon loops giving 



where 



In— + 2F(m/A)J (26) 

with jj, the renormalization scale and 

- j V^ 2 ~ 1 ta nh^ 1 V x 2 - 1 ; x < 1 , , 

\ .^i^Ttan- 1 ;x>l 1 ' 

Analogously, the corrections to can be computed. With the corrections mentioned 
above, the B[B / y coupling is now 

f T =^{^T_-iA+ f"W - 9/2) j-r^K,, A 1/2 ) + -^/K A 1/2 )l (28) 



whereas the B\B^j coupling now is 

" r = 7S fei " ir) + ff {Jr/ im " As/2) + *k 1(m «' M (29) 

The quantities g, g', f and /" are expected to be of order one on dimensional grounds. In 
fact the NRQM predicts g = 1. Unfortunately, the D* lifetime has not been measured yet, 



preventing the extraction of g from D* — > Dir. The experimental upper limit [15] is g < 0.7. 



The model of Ref. |TJ| predicts g = 0.32. Taking g to be somewhere in between we see that 



the magnetic coupling in ( p5j) still remains a quantity of order one or possibly larger. The 
same can be said about the corrections to the contributions from axial-vector mesons, where 
the corrections do not necessarily reduce the value of the couplings given the different signs 
in (|28|) and as well as the potential relative phases between these strong couplings. 
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So far we have not made a distinction between the contributions coming from the B[ 
axial-vector meson, corresponding to ji = 1/2, and the Bi corresponding to ji = 3/2. In 
general one can think that their relative contribution is only governed by the size of 
and an d their respective decay constants. There is already evidence from the charm 

meson system that the (0 + , 1 + ) is very broad relative to the (1 + , 2 + ). This is mostly due to 
the fact that the latter does not couple to the ground state to leading order but only through 
a D-wave amplitude [|1^]. Therefore the narrow width approximation implicitly used in the 
calculation of (|T0|) should be reconsidered. In principle, we expect a modest suppression of 
the B[ contribution. We have also neglected the possibility of (B[,Bi) mixing which could 
be sizable. It has been pointed out in Ref. [12[ that in the D meson system this indeed may 
occur. 

There will be additional corrections coming from terms suppressed by and by the 
chiral symmetry breaking scale A x , which are sub leading terms in the heavy mass and chiral 
expansions respectively. However we believe that one should not expect these or the one 
loop corrections to reduce the value of the couplings significantly, which therefore will remain 
of order one or larger. In the case of /zy, the B*B^f coupling, the one loop corrections are 
always of the oposite sign. However this is not always true for the /ia's, where the corrections 
could even enhance the value of the coupling considerably. 

Perhaps the most important question in terms of the phenomenological impact of these 
decays is the approximation made in Sec. II, namely the neglect of the form-factor suppres- 
sion that affects /iy and \x^_ and that would soften the spectrum. This peaks at a photon 
energy of ~ 1.2GeV when the couplings are approximated to be constant. When an energy 
suppresion is taking into account the average photon energy would drop typically to a few 
hundred MeV. This suppression would result in a smaller value of by a factor that, as 
for the shape of the photon spectrum, strongly depends on the energy dependence chosen. 
Typical energy dependences would reduce R% in ([12]) by a factor of 2 to 4. The first estimate 
is obtained with a monopole type suppression whereas the second case corresponds to an 
exponential suppression. In both cases the energy scale was chosen to be lGeV. 

On the other hand, we lack an understanding of the axial-vector meson decay constants 
and there is a lot of room for theoretical improvements. For the purpose of this letter, 
however, it was sufficient to show that a rough estimate of the plausible values of the 
relevant couplings suggests that the branching ratio for the photonic decays B — > fiwj is 
comparable to B —>■ fiu and it could be even greater. 



IV. CONCLUSIONS 



Radiative leptonic decays of heavy mesons are very interesting, not only as a possible 
background to purely leptonic decays, but also because they yield new information about the 
strong and electromagnetic interactions of heavy hadrons. As expected, they compensate for 
their suppression by a factor of a by avoiding helicity suppression. In the charmed mesons, 
D and D s the effect is small for the muon but sizeable for the electron. Allowing for a typical 
form-factor suppression in /iy and and assuming a large range for the axial- vector decay 
constants we have 

« (1 - 8) x KT 2 /iyGeV 2 (30) 
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which would translate into Br(D — > /Uz/7) ~ 10 -5 . 

The situation is very different for the B meson decays. Since the helicity suppression is 
so large, we benefit much more by avoiding it; we get 

R% » (1 - 4)/£ (l + l -{ ll/2 + 73/2 ) 2 ) GeV 2 (31) 

We have seen in Sec. II that the value of /iy in the B mesons is likely to remain a quantity 
of order one even after considering corrections that might reduce it. This indicates that a 
measurement of fs can only be achieved in the decays to muons if the radiative decays are 
correctly substracted. 

This class of decays clearly deserves further study. In particular it requires knowledge of 
the strong couplings g, g', f and /" mentioned above. In any case it seems reasonable to 
expect values in the range Rf^ = (1 — 20). This would translate into 

Br(B -> luj) « (10~ 7 - 10~ 6 ) (32) 

Therefore the observation of these decays would be of great interest, adding to our under- 
standing of /iy, and the decay constants of orbitally excited heavy mesons. 

There are in principle several other decays of heavy mesons that proceed via the SD 
mechanism. For instance, the decay B s — > [i + fi~y will be enhanced over the helicity sup- 
pressed B s — > ji + n~ by an expression similar to ([H|). In the Standard Model, the latter has 
a branching ratio of ~ 10~ 9 whereas the radiative decays could be an order of magnitude 
greater. In some theories beyond the SM these decays are expected to have larger branching 



ratios. For instance, it is argued in Ref. |L6[ that Extended Technicolor scenarios enhance 



Br(B s — > ) by up to two orders of magnitude, which would imply, together with 
Br(B s — > /i + /i~7) m 10~ 7 — 10~ 6 . This puts this radiative decay within experimental reach 



and, together with B — ► X [17], can impose severe constraints on physics beyond the 
SM. 

On the other hand, B — ► D^ylv and any other exclusive semileptonic decay with a 
photon added to the final state can be considered. Similar results may hold for non-leptonic 
decays. Although these decays are indeed suppressed by a relative to the corresponding 
non-photonic decays, some of them are of interest in their own right and might became 
observable in the near future. 
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FIGURES 



Fig. 1: Bremsstrahlung diagrams entering in B — > jj,i"y. The black squares denote the 
action of the weak current. 

Fig. 2: Structure Dependent diagrams entering in B — ■> /iz/7. The black circles denote 
the action of the couplings /iy and /i^. (a): Vector meson contribution, (b): Axial- Vector 
meson contribution. 

Fig. 3: One Loop corrections to /iy and 1$ . Dashed lines denote goldstone bosons while 
solid lines denote heavy mesons. 
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